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At a generic quantum critical point, the thermal expansion a is more singular than the specific heat 
c p . Consequently, the "Griineisen ratio", F = a/c p , diverges. When scaling applies, T ~ r p~ 1 /( vz ) 
at the critical pressure p = p c , providing a means to measure the scaling dimension of the most 
relevant operator that pressure couples to; in the alternative limit T — > and p ^ p c , Y ~ ^ with a 
prefactor that is, up to the molar volume, a simple universal combination of critical exponents. For a 
magnetic-field driven transition, similar relations hold for the magnetocaloric effect (l/T)dT/dH\ s . 
Finally, we determine the corrections to scaling in a class of metallic quantum critical points. 

PACS numbers: 71.10-Hf, 71.27.+a, 75.20.Hr, 71.28.+d 

Introduction: The anomalous behavior observed in an 
increasing number of systems, ranging from insulating 
magnets and heavy fermion compounds to cuprate su- 
perconductors, has been attributed to the presence of 
quantum critical points (QCPs). These occur in systems 
where a continuous quantum phase transition (QPT) at 
T = is induced by tuning some control parameter 
like pressure p, doping or magnetic field H. Such zero- 
temperature critical points can determine the properties 
of materials in a wide range of temperatures. In gen- 
eral, quantum critical points are more difficult to char- 
acterize compared to their classical counterparts. At a 
classical critical point, thermodynamic quantities typi- 
cally diverge; the associated critical exponents histori- 
cally played a central role in our eventual understanding 
of scaling and universality. Some of these divergences, 
however, have to disappear at a QCP: there are con- 
straints placed by the third law of thermodynamics due 
to the very fact that the transition takes place at zero 
temperature. Here we show that the Griineisen ratio [1,2] 
diverges at any QCP, in a way that provides a novel ther- 
modynamic means of probing quantum phase transitions. 
We define the Griineisen ratio Y [1-3] in terms of the 
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FIG. 1. Schematic phase diagram with a QCP. 
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where S is the entropy and V m — V/N the molar volume. 
In ordinary situations, pressure dependences are regular 
and a finite Griineisen ratio is expected as is indeed ob- 
served in all previous measurements of this quantity in 
the literature. Such a regular dependence is typically de- 
scribed by assuming that the system is dominated by a 
single energy scale E* (e.g. the Fermi energy in a metal 
or the Debye energy if acoustic phonons dominate) so 
the molar entropy takes the form S/N = f(T/E*). The 
Griineisen ratio is then temperature independent and 
given by [1,2,4,5] T — v 1 E , . However, this formula 



already suggests that a diverging Y can be expected when 
some energy scale E* vanishes as it happens at a QCP. 

Divergence of the Griineisen ratio at QCPs: A quan- 
tum critical point is reached in a singular fashion by tun- 
ing some external parameter and, in general, this external 
parameter is thermodynamically coupled to pressure. In 
the low temperature limit, the singular terms of S and 
T in Eq. (1) cancel out leaving T to depend only on 
singularities associated with the pressure p. As the pres- 
sure controls the QPT, such a singularity always exists 
and the Griineisen ratio diverges at any QCP. This diver- 
gence is entirely determined by the scaling dimension of 
the control parameter, which is the most relevant opera- 
tor to which the pressure couples. As shown below, this 
leads to a T dependence, Y ~ l/T 1 /^ [see Eq. (8)]. In 
other words, the temperature exponent of the Griineisen 
ratio provides a direct means to measure vz and, as a 
result, characterize a QCP. Put in a slightly different 
way, the thermal expansion contains valuable information 
complementary to that obtained from the specific heat: 
while c p measures the response to T (y-axis in Fig. 1), 
a describes the response to the tuning parameter of the 
QPT, the second relevant variable at a QCP (x-axis in 
Fig. 1). This has to be contrasted with a classical phase 
transition. There, generically, only one relevant operator 
exists to which both T and p couple. Accordingly, Y will 
be constant close to a classical transition. 
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To observe the singular behavior of T or the thermal 
expansion the pressure has to couple sufficiently strongly 
to the critical dynamics. This is for example the case in 
heavy fermion compounds where the intricate competi- 
tion between magnetic interactions and the Kondo effect 
can be tuned by pressure, doping or magnetic field to 
yield a QPT, typically from a metallic antifcrromagnet 
to a metallic paramagnet. The high sensitivity to pres- 
sure arises from the exponential dependence of the Kondo 
temperature on system parameters. Whether the transi- 
tions in these systems conform to the Gaussian picture 
associated with T = spin-density wave (SDW) transi- 
tions [6,7] or are non-Gaussian as in a locally quantum 
critical point [8] is a question of great current interest. 

If the control parameter of the QPT is not pressure but 
an external magnetic field H, the role of the Griinciscn 
ratio is played by the ratio of the T-derivative of the 
magnetization M (per mole) to the molar specific heat, 
for either fixed pressure or fixed volume, 
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It can be determined directly from the magnetocaloric 
effect - the change of temperature in response to an adi- 
abatic (S = const.) change of H. 

In the following, we go beyond these general consid- 
erations by carrying through a) a more detailed analysis 
based on the assumption of scaling and b) a model study 
on the spin-density-wave (SDW) QCPs in metallic sys- 
tems. The latter is a model system that is above or equal 
to the upper critical dimension, so corrections to scaling 
are important; it is of direct interest in comparing with 
experiments in heavy fermion compounds. 

Scaling analysis: Close to any QCP, the correlation 
length £ diverges as a function of a control parameter r, 
£ - \r\~ v , where e.g. r = (p-p c )/p c orr = (H-H c )/H c . 
Correspondingly, a typical correlation (imaginary) time, 
£ r ~ £ z , diverges as the QCP is approached. The "dy- 
namical critical exponent" z depends on the dynamics of 
the order parameter and relates time and length scales. 

If one assumes that the critical behavior is governed by 
£ and £ T (a more careful discussion of this assumption is 
given below), the critical contribution to the free energy 
per mole, F cr = F — T rcg , can be cast into the following 
standard scaling Ansatz (using hyperscaling) 



1 cr 

~N~ 



-Por 



-po 



v{d+z) 



T 



f , 

1 T r v 

(d+z)/z 

f 



(T/T )V(^) 



(3) 



where po and T are non- universal constants, while f(x) 
and f(x) are universal scaling functions. Obviously, 
f(x — ► 0) « / (0) + xf'(0) + ... is regular as there is no 
phase transition at r = 0, T > (see Fig. 1). The limit 



f(x — ► 0) = /(0) + cx Va+1 describes the low temperature 
behavior of the phases to the left or right side of the QCP 
(in general different for r > and r < 0). Note that the 
exponent yo > has to be positive due to the third law of 
thermodynamics. It characterizes the power-law behav- 
ior of the specific heat c p ~ T Va , e.g. yo = 1 for a Fermi 
liquid, yo = 2 for a d-wave superconductor in d = 2, or 
yo = d and d/2 for an insulating antiferromagnet and 
ferromagnct, respectively. 

Thermodynamical quantities are easily obtained from 
(3). The critical contribution c cr to the specific heat at 
r = is given by 



c cr (T,r = 0) = 
and for T -> 0, r ^ 
c cr (T^0,r) 
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Similarly, in the case of a pressure tuned QCP with r = 
(p — p c ) /p c the critical contribution a cr to the thermal 
expansion reads 



a cr (T,r = 0) = 
and for r ^ 
a cr (T^0) = - 
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The thermal expansion is more singular than the specific 
heat leading to a Griinciscn ratio 



T cr (T,r = 0) 



where the prefactor Gt 
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some non- universal parameters (p c and To). We reach 
the important conclusion that the temperature exponent 
of the Griineisen ratio is equal to l/vz. 

In the other limit T — > 0, r ^ we obtain the universal 
result 



T cr (T^0,r) = -G r 



1 
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Remarkably, even the (generally unknown) scaling func- 
tions cancel out in the amplitude G r , leaving only a com- 
bination of critical exponents and the dimensionality: 
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v(d - ypz) 
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(10) 



Note that the universality of this prefactor is connected 
to the third law of thermodynamics - a finite residual 
entropy per volume (yo = 0) would spoil this result. 
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It is rather difficult to measure thermal expansion in- 
side a pressure cell. However, in many systems doping 
acts like "chemical pressure" . If doping x and pressure 
p can be quantitatively related, p — p c — c(x — x c ), a 
measurement of T for different samples at ambient pres- 
sure can be used to check the prediction (9) quantita- 
tively. For generic tuning parameters, we need to substi- 
tute (dr/dp) for l/p c in Eqs. (6,7) and modify Eqs. (8,9) 
accordingly. 

Similarly, for a QCP tuned by magnetic field [r = (H — 
H c )/H c ] one obtains for the magnetocaloric effect 
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Again, in the T — > limit the prefactor (10) is universal. 
The T-dependence of Tn,cr at r = is also given by (8). 

It is interesting to compare the above with the case of 
a quantum critical line, where the critical behavior is not 
restricted to a single point but to a finite (pressure) inter- 
val. Here, since only marginal and irrelevant operators 
exist for T = 0, T cr can diverge at most logarithmically 



±logT. 



(12) 



Conversely, if T diverges algebraically for T — > a critical 
line scenario can be excluded. 

Applicability of scaling: The applicability of the scal- 
ing results (3-11) depends on a number of assump- 
tions. Most importantly, in an actual experiment not 
T cr = a cr /c cr but r = a/c is measured and sometimes 
non-critical contributions can dominate (for an example 
see below). To verify Eqs. (8) and (9) in a situation where 
c cr is subleading, one would have to subtract carefully the 
non-critical contributions to the specific heat. 

Generally, the scaling Ansatz (3) holds only below the 
upper critical dimension (d + z < 4 within $ 4 theories). 
At the upper critical dimension, logarithmic corrections 
to scaling arise. Above the upper critical dimension, the 
scaling argument can be spoiled by the presence of so 
called "dangerously irrelevant operators" : the free energy 
is a singular function of irrelevant variables. Explicit 
calculations (see below) for the case of an SDW transition 
[6,7] show that on the paramagnetic side the irrelevant 
operator at most leads to logarithmic corrections. 

A more subtle question is whether one of the basic as- 
sumptions underlying the scaling approach (3) holds: Is 
there a single diverging time-scale close to the QCP? For 
example, in a nearly magnetic metal the answer to this 
question is not obvious as there are at least two types 
of low-energy degrees of freedom: magnetic fluctuations 
and fermionic quasiparticles [9-11]. This can indeed lead 
to a breakdown of simple scaling relations as shown e.g. 
by Belitz et al. [9] . In the case of a local critical point in- 
duced by a (non-local) magnetic transition, as has been 
suggested by one of the authors in [8], two scaling di- 
mensions need to be considered: one associated with the 



tuning of the long-wavelength fluctuations and the other 
with the tuning of the local fluctuations. 

SDW transitions: We now turn to more specific calcu- 
lations at SDW quantum critical points, for two reasons. 
First, they allow us to address a number of questions con- 
cerning the scaling results: How do corrections to scaling 
arise at the upper critical dimension? Are the scaling re- 
sults valid above this dimension? What happens, if the 
prefactor d — y^z in Eqs. (9,11) vanishes? Second, our 
calculations are important for the purpose of assessing 
the relevance of SDW QCPs to the magnetic quantum 
phase transitions in heavy fermion compounds. 

Our starting point is the Ginzburg-Landau-Wilson 
functional of Hertz [6] : 

sm = E (*+? 2 + ^W-J 2 + s (4) , (is) 

S (4) = u J o dT J ddr iHr,r)}\ 

with T q = T q z ~ 2 , where z = 2 for an antiferromagnetic 
SDW transition in a metal. The z = 3 theory may be 
used to describe the critical endpoint of a metamagnetic 
first order transition [12] in d = 2, 3. In the case of a 
ferromagnetic QCP in d = 3, the model (13) with z = 3 is 
only valid up to logarithmic corrections and breaks down 
in d = 2 [13]. For commensurate 2D magnetism coupled 
to 2D fcrmions there arc additional singularities in the 
fermion-collective-mode coupling [11]; these singularities 
are absent when the fermions are taken to be 3D [14]. 
Following the renormalization group scheme adopted by 
Millis [7,15] we have calculated the thermal expansion 
and the Gruneisen ratio for d = 2, 3 and z = 2, 3 on the 
non-magnetic side of the phase diagram, S > 6 C . Details 
of the calculation will be reported elsewhere. 
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TABLE I. Resuits for SDW-QCPs in the Fermi iiquid 
regime r = S — S c ^> T 2 ' z . For a pressure tuned QCP, one 
obtains T cr = {dr / dp)T TiCr /V m using r = (p — p c )/Pc, and 
^H,cr = (dr /dH)Y rtCr for r = (H — H c )/H c . Non-universai 
prefactors of a cr and c cr are not shown. The prefactors of 
r cr and Th,ct are (up to the iogarithmic correction for d = z) 
universai. Note that for d = 3, z — 2 the specific heat is 
dominated by a non-criticai contribution c p ~ T. 

The results are summarized in Tables I and II. Up to 
logarithmic corrections the results obey the scaling forms 
(3-11) with v = 1/2, yo = 1. Note that for d = z the pre- 
factor in (7) and (9,11) vanishes. The 1/r dependence of 
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TABLE II. Results for SDW-QCPs in the quantum critical 
regime r = 8 - S c < T 2/z (cf. table I). 



a cr for d = z arises from a T 2 logl/r correction to F cr 
not captured by scaling. For the quantum critical regime 
in d = = 2 the thermal expansion is logarithmic. 
The argument of the logarithm is a power of T for d + 
z > 4 and is itself logarithmically dependent on T for 
d+z = 4; these features reflect the dangerously irrelevant 
or marginal nature of the quartic coupling u. 

In addition to the critical contributions, the measured 
quantities also contain non-critical background compo- 
nents. We list here the full results for the purpose of com- 
parisons with experiments in heavy fermion compounds 
undergoing an antifcrromagnetic transition (z — 2). Con- 
sider first d = 3. At the QCP (r = 0) 



a = ai T 1/2 + a 2 T , 



(14) 



where the a 2 term comes from the (fermionic) back- 
ground contribution. However, approaching the QCP in 
the Fermi-liquid regime 



a={a 1 /r 1 / 2 +a 2 )T. 
For d = 2 and z = 2, we have at the QCP (r = 0) 



a = ailog[&log ^] + ft 2T , 



(15) 
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and in the Fermi-liquid regime approaching the QCP: 



a = (ai/r + a 2 )T . 



(17) 



In two dimensions, the thermal expansion at r = 
diverges in the zero temperature limit in sharp contrast 
to the textbook statement that a(T — > 0) = 0. Still, it is 
straightforward to show that our results satisfy the third 
law of thermodynamics. As a = —(l/V)dS/dp, we can 
write, for generic pressure, 



S(p, T) = S(p c , T) - aV dp - [ P aV dp , 



(18) 



where characterizes the crossover between the QC and 
FL regimes. S(p,T — > 0) — > due to a vanishing inte- 
gration region [(p* — p c ) oc T] over which a is divergent. 

We now briefly discuss the experimental implications 
of our results. Many heavy fermion compounds have long 
been known to show a Griineisen ratio that increases to a 



very large value as temperature is lowered [16,5]. Exper- 
iments are also becoming available in the heavy fermion 
metals tuned to an antiferromagnetic quantum critical 
point, making possible a systematic comparison with our 
theory [17]. It is hoped that the present paper will stim- 
ulate similar measurements in other kinds of (real and 
putative) quantum critical materials. 

In conclusion, we argue that the Griineisen ratio and 
the magnetocaloric effect are divergent at any QCP. In 
addition, they can be used to measure the scaling di- 
mensions and to check the very existence of a quantum 
critical point. 
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